Dynamics of a thin liquid film with surface rigidity and spontaneous curvature 
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The effect of rigid surfaces on the dynamics of thin liquid films which are amenable to the lubri- 
cation approximation is considered. It is shown that the Helfrich energy of the layer gives rise to 
additional terms in the time-evolution equations of the liquid film. The dynamics is found to depend 
on the absolute value of the spontaneous curvature, irrespective of its sign. Due to the additional 
terms, a novel finite wavelength instability of flat rigid interfaces can be observed. Furthermore, the 
dependence of the shape of a droplet on the bending rigidity as well as on the spontaneous curvature 
is discussed. 
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INTRODUCTION 

Surfactant mono- and multilayers at a free interface, 
which are most commonly formed by lipid molecules, ex- 
hibit extremely rich phase diagrams with an abundance 
of different ordering states []], Q • Besides their fascinat- 
ing thermodynamics, they have been intensively stud- 
ied as model systems for biological membranes which 
are composed of phospholipid bilayers Furthermore, 
floating surfactant layers are important for technical ap- 
plications as they can be transferred onto solid substrates 
utilizing, e.g., Langmuir-Blodgett transfer, allowing for a 
controlled coating with an arbitrary number of molecular 
layers and even for the self-organized nano-patterning of 
surfaces 

Depending on the chemical and thermodynamical 
properties of the used material, surfactant mono- and 
multilayers are rigid to some extent. The consequences 
are twofold: On the one hand, the layer has a certain 
bending rigidity, so that additional work has to be sur- 
mounted in order to deform the surface. This bending 
rigidity is measurable by x-ray scattering and depends 
strongly on the surfactant material and the thermody- 
namic phase of the layer [8, 9]. On the other hand, the 
layer can exhibit spontaneous curvature due to different 
chemical properties of the head- and tailgroups of the 
surfactant molecules [Tq| . 

One can expect these effects to be even more pro- 
nounced if the liquid support of such a surfactant layer 
becomes a very thin film, as is the case for example in 
coating processes, so that interfacial influences become 
increasingly important. The dynamics of thin liquid films 
is accurately described within the framework of the lubri- 
cation approximation |ll|, [l2| , which has been used exten- 
sively throughout the literature to describe films covered 
with layers of soluble- and insoluble surfactants 13l-[l5|. 
However, the focus of these investigations has been on the 
effect of surfactant induced surface tension gradients or 
on monolayer thermodynamics, whereas surface rigidity 
is commonly neglected. It is the purpose of the present 
paper to generalize the lubrication approximation of thin 
liquid films in this respect by incorporation of the bend- 




substrate 

FIG. 1. Schematic drawing of a thin liquid film with a rigid 
surface layer. The shape of the film is described by its height 
profile h(x,t), which gives the film thickness above point x 
at time t. 



ing rigidity and the spontaneous curvature of a surface 
layer. Following this approach, we shall find that the 
surface rigidity gives rise to a pair of antagonistic effects: 
While short wavelength fluctuations are suppressed by 
the bending rigidity, the spontaneous curvature can lead 
to instabilities at intermediate wavelengths. 



THIN LIQUID FILMS WITH RIGID SURFACE 
LAYERS 

In the spirit of Helfrich, the energy of a rigid interfa- 
cial layer can be written as the sum of integrals over the 
mean curvature and the Gaussian curvature of the sur- 
face |To|, [l6|. According to the Gauss-Bonnet theorem, 
the latter integral is a topological quantity, dependent 
only on the genus p of the considered surface S. De- 
noting the principal curvatures of S by K\ and k<i, the 
spontaneous curvature by kq, and the surface tension by 
a, the energy E Y due to the rigidity of a surface S is 
accordingly obtained as 



E r = ^ JdS (*i 



+ «2-«o) 2 + 27r*; g (l-jp). (1) 



The positive material constants k c and k g are the bend- 
ing rigidity and the modulus of the Gaussian curvature, 
respectively. Here, we will focus on thin films which are 
entirely covered by a homogeneous surfactant layer with- 
out defects, that is, the case of constant surface topology. 
If the liquid film is described by its height profile h(x,i) 
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as shown in figure [TJ the Helfrich energy can be written 
a a functional E r = E r [h] . 

In equilibrium there has to be a balance of forces SW = 
at the surface, where 



SW = JdSSh jp-a(^+^ 2 )-^ d + ^J. 



(2) 



Here p is the pressure in the liquid phase and 7Td denotes 
the disjoining pressure describing the interaction between 
substrate and liquid which has to be taken into account 
for very thin films. To calculate the functional derivative 
of E T , we use the fact, that one can rewrite the mean 
curvature as k\ + ^2 = —V • fi. where 
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with a 



(V/1) 2 



is the unit normal of the water surface, pointing into the 
vapor phase. We find 
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(3) 



Introducing the scales ho, Zo, and to for height, length, 
and time, one can derive the time evolution equation 
for the nondimensionalized height profile H(X,T) := 
h(x/lo,t/to)/ho within the lubrication approximation. 
To this end, we follow precisely the steps in chapter II. B 
of the derivation by Oron et al. with the difference, 
that our force balance equation contains the additional 
contribution 5E r /5h. Expanding the functional deriva- 
tive (j3j) in powers of the parameter e = ho/lo, which is 
small for thin liquid films, one obtains 



SE T k c . { 2e . enl ^ TT 

ih=^-U VAH -^ H 



4e 2 ^ Q 
I 2 



(AH)VH - -V(ViT) 2 



- [HAH)(VHf 
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Inserting this result into the force balance equation (|2]) 
and expanding the remaining terms in powers of e as well, 
we obtain, after dropping terms of higher orders in e, the 
following formula, which is a generalization of equation 
(2.24b) in [Tl|: 



-p + n = C^AH 



C^A'H, 



where P = ehop/d^Uo) and II = eho^d/ {l^Uo) with char- 
acteristic velocity Uo = lo/to are the nondimensionalized 
pressure and disjoining pressure and 



0- = 



UofJL 



and C T 



UofJL 2/ 2 



(4) 



With this result at hand, one can follow the remaining 
steps of the standard derivation in to obtain the time 
evolution equation 

d T H = -V • j-^V [C^AH - C^A 2 H - 11(H)] J . 

(5) 

It is a noteworthy result, that the thin film dynamics does 
not depend on the sign of the spontaneous curvature, but 
only on its absolute value. 

Typical values for the bending modulus and the spon- 
taneous curvature are k c w 10 -19 J [1, 0] and \k>o\ ~ 
10 7 m _1 - 10 9 m _1 respectively. One can thus 

roughly estimate k c n^/2 w 10 _5 J/m 2 - 10 _1 J/m 2 . 
Comparing this range of values to the surface tension of 
water in absence of any surfactant, <r a b s ~ 72- 10 -2 J/m 2 , 
we note that the contribution due to spontaneous cur- 
vature can be of the same order of magnitude or even 
larger. The latter implies that spontaneous curvature 
can lead to a negative effective surface tension. As- 
suming Iq « 10 _T m - 10 -6 m one can further estimate 
k c /(2ll) w 10- 7 J/m 2 - 10- 5 J/m 2 . Thus, Cf/ will in 
most cases be orders of magnitude smaller than Cf 1 . 



LINEAR STABILITY OF FLAT FILMS 

Any homogeneous film profile H(X,T) = H = const, 
solves equation (j5j). Considering small perturbations 
r)(X, T) of a flat film, we insert H(X, T) = H + rj(X, T) 
into equation (j5j) and keep only terms linear in 77. Using 
the plane wave ansatz 77 ~ exp(ifc-X + AT) the dispersion 
relation A(fc) is obtained as 



X(k) = -^ 



fl' + Cf 1 *; 2 



Cf/fc 4 



(6) 



where we use the shorthand notation II 7 := 
dH~R(H)\ H _£ and k := \k\. This result complies 
with the physical intuition, that, for vanishing sponta- 
neous curvature, the bending rigidity will damp short 
wavelength fluctuations and thus has a stabilizing 
influence. However, if C^ 1 takes negative values, this 
can lead to an instability for intermediate wavenumbers: 
For Cf 1 < and tl f = Cf 2 /(4Cf I 1 ) the eigenvalue 

A(fc cr ) of the mode k cr = 



-Cf 1 /(2Cf I 1 ) crosses 
the imaginary axis (solid line in figure [2]), so that for 
< n / < Cf 2 /(4C n 1 ) and Cf 1 < there is a band of 
unstable modes k\ < k < k<i (dash-dotted line in figure 
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FIG. 2. Dispersion relation X(k) (see equation ©) for four 
different parameter sets: Cf 1 = -0.1, Cf 1 = 0.1, ft' = 0.025 
(solid line), Cf 1 = -0.1, Cf 1 = 0.1, ft' = 0.0148 (dash- 
dotted line), Cf 1 = 0.1, Cf 1 = 0.1, ft' = -0.0157 (dashed 
line), and Cf 1 = 0.1, Cf 1 = 0.1, ft' = 0.0148 (dotted line). 




FIG. 3. Snapshot of a numerical simulation of the breakup of 
a homogeneous film with H — 1.2, Cf 1 — —0.1, Cf 1 — 0.1. 
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Obviously k\ = for negative Ir and arbitrary values 
of Cf 1 (dashed line in figure [2]). The wavenumber & max 
corresponding to the maximum of A(fc), that is, the most 
unstable mode, is given by 



\ 



H \- c T +V C i " 3C n n ' 



81C n 2 



6Cf I 1 n / - Cf 2 + CfVCf 2 - 3Cf I 1 n / . 



with 



A(&max) 



For f[ > 0, Cf 1 > 0, and Cf 1 > 0, the whole spectrum 
is damped (dotted line in figure [2]). 

A numerical simulation [18| of a randomly perturbed 
fiat film of height H = 1.2 with Cf 1 = -0.1 and 
Cf 1 =0.1 was carried out using the disjoining pres- 
sure U(H) = AH~ 3 (l - H~ 3 ) with Hamaker constant 
A = 6.5 • 10 -2 . This is precisely the situation corre- 
sponding to the dash-dotted line in figure [2j The fiat 
film quickly breaks up into small droplets which then un- 
dergo a process of coarsening, as is commonly observed 
after film rupture. 

Remarkably, however, these results predict, that rigid 
layers with spontaneous curvature can be used to destabi- 
lize very thin films, which would remain flat in absence of 
the surfactant layer, and to obtain droplets of extremely 
small volumes after film breakup. 
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FIG. 4. The different shapes of a droplet obtained for Cjj 1 = 
0.1 and Cf 1 = (dash-dotted, line), Cf 1 = -0.05 (solid line), 



Cf 



-0.07 (dotted line), Cf 



-0.1 (dashed line). The 



overall volume of liquid was equal in all four cases. The image 
shows a closeup of a steady state on a periodic domain of size 
L = 100. 



THE SHAPE OF A DROP COVERED WITH A 
RIGID LAYER 



It was shown in the preceding section, that surface 
rigidity affects the linear stability of flat films. Now, we 
want to investigate how the rigidity affects the shape of 
the droplets that are formed during the rupture of an un- 
stable film. To this end, we have considered steady state 
solutions of equation || on a one-dimensional domain 
which were obtained by relaxation, that is, direct numer- 
ical integration of equation (0 . When a steady state was 
obtained for one set of parameters, it has been used as 
initial condition for the subsequent run. Thus, the over- 
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droplet obtained for C^ 1 = 
10 (dash-dotted line), 
jj — uu ^ciDiic»a imc; , v^jj — 50 (solid line). The overall 
volume of liquid was equal in all four cases. The image shows a 
closeup of a steady state on a periodic domain of size L = 100. 



FIG. 5. The different shapes of 
0.1 and C" 1 = (dotted line), 
CVt 1 = 30 (dashed line), C^ 1 = 



all volume of liquid within the integration domain was 
constant throughout a whole parameter scan. 

First, Cfj 1 = 0.1 was held constant while Cj was var- 
ied from to —0.1. The obtained shape of the droplet 
depends sensitively on Cj - , as is shown in figure 2J The 
drop becomes steeper and develops more and more pro- 
nounced undershoots below the precursor height H = 1 
as the parameter is decreased. 

In a second parameter scan Cf 1 = 0.1 was fixed while 
C^ 1 was varied from to 50, as is shown in figure [5] The 
results indicate that the scaled bending modulus k c / (2/q) 
needs to be 10 - 100 times higher than the effective sur- 
face tension in order to change the shape of the droplet 
significantly. 



CONCLUSIONS & OUTLOOK 

A generalization of the equations of thin film flow to 
the case of rigid surface layers with spontaneous curva- 
ture has been presented. Our results show, that the spon- 
taneous curvature can lead to a significantly lower and 
even negative effective surface tension, yielding a new 
finite wavelength instability resulting in the breakup of 
very thin flat films. It turned out, that only the absolute 
value and not the sign of the spontaneous curvature is 
important in the lubrication regime. 

We have shown that the sensitivity of steady droplet 
solutions on C^ 1 is very low. Combined with the fact 
that it is very small in most physical situations, we con- 
clude that its most important contribution is the stabi- 
lization of short wavelength fluctuations in the case of 



negative Cf . 

It would be of interest to further extend the model to 
the case of surfactant layers with evolving topology, that 
is, to layers that can break up and form holes, leading 
to contributions <~ k g to the surface energy. This could 
also be applicable to binary surface layers consisting of 
material in two different thermodynamical phases, if only 
one phase has significant rigidity. 

This work was supported by the Deutsche Forschungs- 
gemeinschaft within SRF TRR 61. 
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